Revisit of Coulomb effects on tt /tt^ ratio in heavy ion collisions 
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We show that the large tt^/ti^ ratio at low pion energies recently observed in 
the reaction of Au+Au at Ebeam/A = 1.0 GeV at SIS/GSI is due to the strong 
Coulomb field of the reaction system. 
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It was found recently by the KaoS and FOPI collaborations at SIS/GSI that the -n~ j-n^ 
ratio increases as the pion kinetic energy or transverse mass decreases The ratio is 

found to depend strongly on the mass of the reaction system. For example, in the reaction of 
Au+Au at Ebeam/^ = 1.0 GcV the ratio increases continuously from about 1 at E^^^ ^ 400 
MeV to about 2 at ^ 100 MeV. While in the reaction of Ni+Ni at E^^am/A = 1.8 GeV 
the ratio goes upto only about 1.15 p. More recently, the E866 collaboration also studied 
the tt^/tt^ ratio as a function of the pion transverse energy Et for Au+Au, Si+Au, Si+Cu 
and Si+Al reactions at AGS/Brookhaven [^]. For the first two reactions the tt^/ti^ ratio 
rises to about 2 at low Et, while for the last two systems the ratio is only about 1 in the 
whole energy range. 

The phenomenon has renewed some interests on the study of /n^ ratio in relativistic 
heavy ion collisions. The main question of interest is whether it is a manifestation of the 
well known Coulomb effect or an indication of some new physics. Moreover, it seems also 
necessary to well understand the isospin and Coulomb effects in order to study the more 
interesting phenomena found in the pion spectra of these experiments, such as, the enhance- 
ment of low/high energy pions. For this purpose we perform a quantitative study on the 
7r^/7r+ ratio in heavy ion collisions at SIS/GSI energies using a hadronic transport model 
|0J^. We show that the large tt^/ti^ ratio found at SIS/GSI is indeed a manifestation of 
the Coulomb effect. Although we consider here reactions at SIS/GSI energies only, results 
of this study may also shed some light on understanding the similar effect found at AGS 
energies. 

For completeness we first mention some of the early studies on the tt^/tt^ ratio in nuclear 
reactions. A large tt^/ti^ ratio at very low laboratory energies was observed a long time 
ago in experiments using protons, alphas and cosmic rays, and was explained as due to the 
Coulomb interactions. These studies were summarized and discussed in great detail in ref. 
0. Later on, much interest was generated by the discovery of structures in charged pion 
spectra at the Bevalac. A very distinctive feature observed in collisions of nearly equal mass 
projectile and target nuclei is an enhancement of the vr"*" yield in the mid-rapidity region, 
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which appears as a peak in the vr"*" spectra at 9cm = 90°. Such a peak was first observed 
at transverse momentum pt ~ 0.4m^c in the reaction of Ar+Ca at Eheam/^ = 1-05 GeV by 



Wolf et al. and at a shghtly higher pt by Chiba et al. |Tl| in the reaction of Ne+NaF at 
Ebeam/A = 0.8 GcV. Additional structure in the low energy region at ^ = 0*^ in the projectile 



frame was found first by Benenson et al. []T2| in the reaction of Ne+NaF at Ef,eam/^ < 400 
MeV. Corresponding to the strong vr^ peak at pt = a. dip was seen in the 7r+ spectrum in 
the same experiments. Some systematic studies on the tt'/tt^ ratio were also performed at 
the Bevalac 1T|-|16 



On the theoretical side, all studies seem to indicate the Coulomb nature of the phe- 
nomenon. By comparing the 7r+, tt^ wave functions in the electromagnetic field generated 
by the moving projectile fragment, Bertsch was able to explain the peak in the tt^/ti^ ratio 
found by Benenson et al. and its beam energy dependence [|12],|l7|. Libbrecht and Koonin 



showed that the tt'^ peak at 6cm = 90° and pt = 0.4m^c could be a result of the Coulomb 



focusing due to the relatively moving charged fragments |[T8|. In their analysis the charge 
distribution was approximated by spreading part of the protons along a line between the 
nuclear centers with the remainder continuing with the unaltered velocities. The classi- 
cal equations of motion for pions were then sloved in the Coulomb field. Gyulassy and 
Kauffmann derived various approximate expressions for ti^/ti^ ratio in the Coulomb field 
of projectile and target remanants as well as the thermally expanding fireball [1^. Later, 
the 7r^/7r+ ratio was also studied in ref. using a statistical model for pion production at 



subthreshold energies 

In the present study we use the hadronic transport model for relativistic heavy ion 
collisions. The model uses explicit isospin degrees of freedom for all hadrons involved in 
the reaction. In particular, isospin-dependent hadron-hadron collision cross sections and 
branching ratios for the decay of baryon resonances are fully incorporated in the model. 
The time dependent, selfconsist mean field for baryons including the Coulomb field and 
an isosymmetry term is a basic ingredient of the model. It has been rather successful 
in studying many aspects of relativistic heavy ion collisions. We refer the reader to our 



previous publications for more details of the model [^. In view of the fact that in most of 
the previous theoretical studies on the subject one had to assume an idealized, simple and 
often static charge distribution during the reaction, the present study has the advantange 
of incorporating dynamically effects of the varying charge distribution. 

Fig. 1 shows the calculated 7r^/7r+ ratio as a function of pion kinetic energy with and 
without Coulomb interactions for the reaction of Au+Au at Ebeam/^ =1.0 GeV. In the 
case without Coulomb interaction the ratio is a constant of about 1.9 for Ekin < 500 MeV, 
beyond 500 MeV the ratio has large statistical error bars. The constant tt^/tt^ ratio is very 



close to the estimate using the isospin argument |23 



I_^5JP + NZ 
' ~ 1+ 5^2 + NZ ' 

which gives 1.95 for the Au+Au reaction. While in the case with Coulomb interactions the 
ratio is larger than the isospin average for Ekin < 150 MeV and decreases as the pion energy 
increases. The cross over of the two calculations at Ekin ~ 150 MeV is a direct result of the 
competition between the Coulomb momentum impulse and the phase space distortion. This 
can be understood qualitatively using classical considerations similar to that of refs. ||T^,|T^ 



In terms of the inclusive cross section cr^ = d^a^/dp^ in the absence of the Coulomb field 
and the change of the phase space density d'^po/d^p, the cross section in the presence of the 
Coulomb field cr^ = d^a^ /dp^ can be written as 



where p is the momentum of pions after moving through the Coulomb filed starting with 
an initial momentum po. In the last approximation we factorized out the isospin factor 
since the /n^ ratio is almost a constant in the absence of the Coulomb interactions. The 
Coulomb impulse is 

6p = ±{p - po{p)) (3) 

for 7r+ and vr^, respectively. Thus, the impulse results in an increase (decrease) for 7r+ (tt^) 
in the (Tq term of eq. (H) which has approximately an exponential form exp{—pQ/2m.^T), 



where T is the slope parameter of the pion spectra. The phase space distortion \d^po/d^p\ 
generally results in a decrease (increase) for 7r+ (vr^), its exact form depends on the dynamics 
of the nuclear charge distribution and the trajectory of individual pions. An example of pion 
emission from the origin of a sphere with the total charge Ze and a radius r is instructive, 
though not very realistic. The potential energy of ± charged pions at the origin is Vq = 
±|Za/r, energy conservation leads to 



p,{p)=pilTpl/p^)y\ (4) 



where pc = \/2m\Vo\. One therefore obtains that 

(IT pI/p'Y^', (5) 



and 



Q3p 



a^{p) = /± ■ aoip) ■ exp{±\Vo\/T)ilTpi/pT- (6) 
Keeping the first order in Za, the /n^ ratio is approximately 

The competition of the two factors m/p'^ and 1/T causes the ti~ /tt^ ratio to change from 
R > Rj to R < Rj as the energy increases to above E^in = T/2. The variation of the -n^ j-n^ 
ratio in the transport model calculations is therefore understandable. 

We now perform a comparison with the experimental data. Before comparing the ratio 
we first compare in Fig. 2 and Fig. 3 the vr"*" and vr" spectrum, respectively. The data were 
taken for the Au+Au reaction at Qi^h = 44'^ by the KaoS collaboration at SIS/GSI ||2|. It 
is seen that the model can very well reproduce the available experimental spectra. Fig. 4 
shows a comparison of the tt^/tt^ ratio. The model calculation is in good agreement with 
the available data at low energies, more complete data will be available from both the KaoS 
and FOPI collaborations at SIS/GSI. It is interesting to mention that a model calculation 
for the reaction of Ni+Ni at E/A = 1.8 GeV, which is almost isosymmetric, shows that the 



'K~ I'K'^ ratio is about constantly 1 within statistical error bars. Comparing to the Au+Au 
reaction, it further indicates the strong Coulomb effect on the tt^/tt"^ ratio in heavy systems. 

In summary, using a hadronic transport model we have shown that the recently observed 
large t:~ /n'^ ratio at low energies in relativistic heavy ion collisions is a manifestation of the 
strong Coulomb field of the reaction system. 
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FIGURE CAPTIONS 

Fig. 1 The calculated 7r~/7r+ ratio with and without Coulomb interactions for the 
reaction of Au+Au at £^(,eam/^ = 1.0 GeV. 

Fig. 2 Comparison to the tt^ spectrum measured at Oiab = 44° ± 4° in the reaction of 
Au+Au at -EfeeamM = 1-0 GcV 

Fig. 3 Same as in Fig. 2 but for the preliminary tt" spectrum. 

Fig. 4 The 7r~/7r+ ratio for the reaction of Au+Au at £^(,eam/^ = 1-0 GeV. The preliminary 
data was measured at Oiab — 44° ± 4° and the calculation is done at Oiab — 45° ±15°. 
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